function test4( )
%TEST4 Summary of this function goes here
%   Detailed explanation goes here

    A = imread('../TestImg/storm.png');
    A = double(A(:,:,1));
    
    disp('A-Size');
    disp(size(A));
    
    figure(1);
    subplot(2,2,1);
    imshow(A./255);
    
    
    psf  = fspecial('laplacian');
    disp('PSF-Size');
    disp(size(psf));
    figure(1);
    subplot(2,2,2);
    imshow(psf);
    
    H = conv2(A,psf,'full');
    disp('H-Size');
    disp(size(H));
    disp('H');
    disp(min(H(:)));
    disp(max(H(:)));
    
    HShow = (H-min(H(:)));    
    HShow = HShow./max(HShow(:));
    
    disp('HShow');
    disp(min(HShow(:)));
    disp(max(HShow(:)));
    
    figure(1);
    subplot(2,2,3);  
    imshow(H);
    
    HRec = MyDeconv2D_2(H,psf, size(A,1)-4, size(A,2)-4, 100);
    
    HRecShow = (HRec-min(HRec(:)));
    HRecShow = HRecShow./max(HRecShow(:));
    disp('HRec');
    disp(min(HRec(:)));
    disp(max(HRec(:)));
    
    figure(1);
    subplot(2,2,4);
    imshow(HRecShow);
    
%     imshow((H_Rec-min(H_Rec(:)))./max(H_Rec(:)));
end

